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ABSTRACT 

Multiple reflection of X-rays may be important when an accretion disc and its hot corona 
have a complicated geometry, or if returning radiation due to gravitational light bending is 
important, or in emission from a funnel such as proposed in some gamma-ray burst models. 
We simulate the effects of multiple reflection by modifying the boundary condition for an 
X-ray illuminated slab. Multiple reflection makes the soft X-ray spectrum steeper (softer) and 
strengthens broad emission and absorption features, especially the K-shell features of iron. 
This may be important in explaining the spectra of sources such as the Narrow-Line Seyfert 1 
galaxy 1H 0707^195. 
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1 INTRODUCTION 

The hard X-ray emission from luminous Seyfert 1 galaxies is of- 
ten modelled as due to an optically-thin corona above a dense ac- 
cretion disc. This model successfully accounts for the reflection 
component often found in the spectra of these objects (Pounds 
et al. 1990). The dense optically-thick accretion disc is often as- 
sumed to be relatively cold, with iron only weakly ionized (e.g. 
George & Fabian 1991; Matt, Perola & Piro 1991; Magdziarz & 
Zdziarski 1995). Sometimes however it may be highly ionized by 
the irradiating X-rays, resulting in a more complex spectrum (Ross 
& Fabian 1993; Zycki et al. 1994; Ross, Fabian & Young 1999; 
Nayakshin, Kazanas & Kallman 2000; Ballantyne, Ross & Fabian 
2001; Rozanska et al. 2002). Generally the disc is assumed to be 
flat. 

Recent numerical simulations of luminous accretion discs 
(Turner, Stone & Sano 2002), particularly when radiation pressure 
is important, suggest that they are clumpy, irregular and possibly 
corrugated. In this case, if the corona lies close to the disc, the solid 
angle subtended by the disc at the corona may be much higher than 
271", the value for a flat disc, and the reflection spectrum itself may 
be reflected. This will certainly happen if the disc is corrugated 
and the corona lies close to the bottom of corrugations. The ob- 
server may then not see any direct emission from the corona but 
only multiply-scattered flux. (Rapid variability of an AGN requires 
that any such corrugations be on a relatively small spatial scale.) 

We attempt here to simulate multiply-scattered reflection 
spectra for ionized discs in a simple generic manner using the code 
described by Ross & Fabian (1993). Multiply-scattered reflection 
spectra have been invoked as one explanation (Fabian et al. 2002) 
for the remarkable spectrum of the Narrow Line Seyfert 1 galaxy 
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1H 0707^495 observed by XMM-Newton (Boiler et al. 2002). They 
might also be relevant if the accretion disc extends well within 6 
Schwarzschild radii so that returning radiation due to light bending 
by the strong gravity of the black hole is important (Cunningham 
1975; Martocchia, Matt & Karas 2002). Multiple reflection is ex- 
pected naturally if the accretion flow consists of clouds rather than a 
disc (e.g. Celotti, Fabian & Rees 1992; Collin-Souffrin et al. 1996; 
Malzac 2001). Multiple reflection could also be important in hy- 
pernova models for Gamma-ray bursts (Meszaros 2001), in which 
the emission is surrounded by ejecta. 

We expect that multiple reflection will enhance both broad ab- 
sorption and emission features in the spectrum, relative to single re- 
flection spectra. The net effect is to produce a considerably steeper 
(softer) spectrum with pronounced broad spectral features around 
the K-shell emission of iron and other elements. 



2 THE MODEL AND METHOD 

As an idealized model, let us consider that the hard continuum orig- 
inates from close to the bottom of a well or hole. If the opening at 
the top subtends an effective solid angle £1 as seen from the bot- 
tom, then a fraction e = Q./2-k of the emitted radiation escapes 
directly; the rest is scattered, absorbed and reflected by the mate- 
rial in the sides and bottom of the well. To a crude approximation, 
this reflected flux has a similar probability of escaping, with the 
rest restriking the surfaces within the well. We model the multiple 
reflection process by simply introducing a partially reflecting mir- 
ror above a flat slab. A detailed model would involve knowledge 
of the precise geometry envisaged, but we expect that our simple 
model will give us a good qualitative guide to what the observer 
would expect to detect. If e is small, an observer may see only the 
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Figure 1. Observed spectra resulting from multiple reflection with escape 
fraction e = 0.25 (solid lines) and from a single reflection (dotted lines) 
for four different values of the ionization parameter §. Spectra for different 
values of § have been offset vertically for clarity. 




Eq. (2) 
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Figure 2. Spectra for £ = 1000 renormalized by dividing the escaping 
flux Fe(E) by the corresponding value at E = 8.8 keV. The solid line is 
the result for multiple reflection with e = 0.25, while the dotted line is the 
result for a single reflection (s = 1). The dashed line shows the approximate 
shape of the Fe K-absorption feature for e = 0.25 predicted using Eq. (2). 



multiply-reflected flux, and it is this reprocessed radiation only that 
we discuss. 

In detail, we treat a surface layer of gas with hydrogen number 
density iih = 10 15 cm -3 and Thomson depth tt = 10. The outer 
surface is illuminated by "direct" radiation in the form of a power 
law with photon index T = 2 and with total flux F (from 1 eV to 
100 keV) corresponding to a specified ionization parameter 



TIH 



(1) 



The penetration of the illuminating radiation and the Compton- 
thick transfer of the diffuse radiation produced via scattering or 
emission within the gas are treated as described by Ross & Fabian 
(1993). At each step in the radiative transfer calculation, the diffuse 
radiation emerging from the outer surface is calculated as described 
by Foster, Fabian & Ross (1986). Only a fraction e of this emergent 
flux, however, is allowed to escape to distant regions. The remain- 
ing fraction (1 — e) of the emergent flux is assumed to be reflected 
back onto the surface, so it is added to the "direct" radiation that 
illuminates the surface. The transfer calculation proceeds until the 
radiation field relaxes to a steady state, at which point the portion 
of the total emergent flux that actually escapes equals the total "di- 
rect" illumination. 



3 RESULTS 

We have calculated multiple-reflection models with direct illumi- 
nation corresponding to £ = 250, 500, 1000 and 2000. For each 
value of the ionization parameter, models were produced with the 
emergent radiation having escape fractions of e = 0.25, 0.50, 0.75 
and 1 .00. Here e = 1 corresponds to escape of all emerging radia- 
tion with no reflection back onto the surface, and these results are 
similar to those of Ross, Fabian & Young (1999). 

Figure 1 shows the observed spectrum when only 25% of the 
emergent radiation actually escapes (and 75% is reflected back) in 
comparison with the spectrum without multiple reflection. Multi- 
ple reflection reduces the escaping flux for E > 1 keV by increas- 
ing the probability that such photons will be absorbed before ulti- 
mately escaping. In addition, it changes the shape of the observed 
spectrum. To make this clearer, Figure 2 shows the same two spec- 



tra for £ = 1000, but with each spectrum renormalized relative to 
the corresponding escaping flux at 8.8 keV (just below the K-edge 
of Fe XXV). The multiple-reflection spectrum exhibits a relatively 
deeper Fe K-edge, a stronger and broader Comptonized Fe Ka line, 
and a steeper spectrum in soft X-rays (~ 1 keV). 

An approximate test can be made of the accuracy of our 
method. Consider a photon energy for which the illuminated ma- 
terial is essentially a pure absorber and coherent scatterer — that is, 
there is no emission due to absorbed or scattered radiation at other 
photon energies. Let A be the albedo for a single reflection by the 
material. Then, if the probability of escape after each reflection is 
e, the net albedo produced by multiple reflection is 

A nct = eA + e(l-e)A 2 + e{l - ef A 3 + ■•• 
eA 



l-d-e)A (2) 

Fig. 2 shows the results of applying this formula to energies 
around the K-edge of iron, where it is approximately valid, for 
e — 0.25. Eq. (2) predicts a slightly deeper absorption feature 
than our multiple-reflection calculation because it does not take into 
account incoherent Compton scattering of higher-energy photons, 
which partially fills in the absorption feature (see Ross, Weaver & 
McCray 1978). 

Figures 3 and 4 provide another way of visualizing the 
changes that multiple reflection produce in the observed spectrum. 
In each case, the ratio of the flux escaping after multiple reflec- 
tions to the flux escaping after a single reflection is shown as a 
function of photon energy. On average, decreasing the escape frac- 
tion e decreases the entire spectrum above ~ 1 keV. As e de- 
creases, however, the broad Fe K-absorption feature decreases by 
a larger amount than the X-ray spectrum as a whole, so that the 
relative strength of this feature is enhanced. On the other hand, 
as e decreases, the broad, Comptonized Fe Ka emission feature is 
suppressed less than the X-ray spectrum as a whole, so that this 
emission feature is also enhanced. At high ionization parameters 
(£ = 2000 and 1000; Fig. 4), multiple reflection suppresses the 
narrow central core of the Fe Ka line, since there is a deep, highly- 
ionized surface layer that Compton scatters line photons that reen- 
ter it. At the lowest ionization parameter (£ = 250; Fig. 3), a de- 
crease in e enhances the Fe XXV line core at 6.7 keV and suppresses 
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Figure 3. Ratio of the observed flux after multiple reflection with escape 
fraction e to the observed flux after a single reflection. The upper panel 
shows the results for £ = 250; the lower panel shows the results for § = 
500. 



the "neutral" iron line at 6.4 keV, since multiple reflection increases 
the effective ionization potential by increasing the total flux inci- 
dent on the surface. In the soft X-ray region, the overall slope of 
the ratio presented in Fig. 3 or 4 steepens as e decreases, indicating 
that multiple reflection steepens the observed spectrum compared 
to that resulting from a single reflection. 



4 DISCUSSION 

Multiple X-ray reflection creates a steep soft X-ray spectrum and 
a strong, broad, iron absorption/emission feature. The power lost 
from the harder energies above ~ 1 keV emerges at lower energies, 
principally in the 0.05-1 keV band. Overall the spectra resemble 
those from slightly higher ionization parameters and much higher 
abundances (Ballantyne, Fabian & Ross 2002). 

A relativistically blurred version of the spectrum for e = 0.25 
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Figure 4. Flux ratios as in Fig. 3, except for £ = 1000 (upper panel) and 
for £ = 2000 (lower panel). 



and £ = 2000 is shown in Figure 5. That has been used to simulate 
what would be observed with the XMM-Newton EPIC pn CCD, 
and this is also displayed in Fig. 5 in the form of a ratio to an in- 
cident power-law of photon index 2.7. The result compares well to 
the spectrum of 1H 0707-495 found by Boiler et al. (2002). If we 
fit the simulated spectrum from 2-10 keV with a power-law plus 
a broad gaussian line, then we find that the best-fitting line has a 
central energy of 4.8 keV, a width of 1.5 keV and an equivalent 
width of about 6 keV. Positive residuals remain over the 6-7 keV 
band. It is not our intention to carry out a detailed fit to the data of 
1H 0707-495, which would require a comprehensive grid of mod- 
els, but merely to confirm that multiple reflection may be relevant 
there, as suggested by Fabian et al. (2002). Since the spectrum of 
1H 0707-495 is unusual, multiple reflection may occur only occa- 
sionally among Narrow Line Seyfert 1 galaxies, and/or the primary 
continuum may be seen in addition to the reprocessed component 
most of the time. 

Multiple reflection can be important when light bending is 
strong enough to cause part of the reflected spectrum to return to 
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Figure 5. Top: Model spectrum for e = 0.25 & £ = 2000 relativisti- 
cally blurred by a kernel appropriate (Laor 1991) for a disk inclined at 
30 deg to the line of sight and extending from 3 to 100 gravitational radii 
(3-100xGM/c 2 ). Bottom: Ratio of fake data obtained by simulating the 
above model for the EPIC pn CCD for 40 ks (and intensity appropriate for 
1H 0707—495) to a a power-law model of photon index 2.7. 



the disk. This occurs close to the centre of disks around maximally 
spinning Kerr black holes (Martocchia & Matt 1996; Dabrowski et 
al 1997; Dabrowski & Lasenby 2001; Martocchia, Matt & Karas 
2002). Strong X-ray reflection has also been invoked (Vietri et al 
2001; Ballantyne & Ramirez-Ruiz 2001) for the production of X- 
ray lines seen in some gamma-ray bursts (e.g. Piro et al 2000). If 
the emission is from a funnel geometry (Rees & Meszaros 2000; 
McLaughlin et al 2002) then the calculations presented here could 
be relevant. Lower abundances are implied than obtained from sin- 
gle reflection calculations. 
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